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ABSTRACT 

Aims. We provide a library of some 7000 SEDs for the nuclei of starburst and ultra luminous galaxies. Its purpose is to quickly obtain estimates 
of the basic parameters, such as luminosity, size and dust or gas mass and to predict the flux at yet unobserved wavelengths. The procedure is 
simple and consists of finding an element in the library that matches the observations. The objects may be in the local universe or at high z- 
Methods. We calculate the radiative transfer in spherical symmetry for a stellar cluster permeated by an interstellar medium with standard 
(Milky Way) dust properties. The cluster contains two stellar populations: old bulge stars and OB stars. Because the latter are young, a certain 
fraction of them will be embedded in compact clouds which constitute hot spots that determine the MIR fluxes. 

Results. We present SEDs for a broad range of luminosities, sizes and obscurations. We argue that the assumption of spherical symmetry 
and the neglect of dumpiness of the medium are not severe shortcomings for computing the dust emission. The validity of the approach is 
demonstrated by matching the SED of seven of the best studied galaxies, including M82 and Arp220, by library elements. In all cases, one 
finds an element which fits the observed SED very well, and the parameters defining the element are in full accord with what is known about 
the galaxy from detailed studies. We also compare our method of computing SEDs with other techniques described in the literature. 
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1. Introduction 

By definition, the rapid conversion of a large amount of gas 
into predominantly massive (> 8Mq) stars, or the result of such 
a conversion, is called a starburst. Starburst galaxies consti- 
tute a unique class of extragalactic objects. The phenomenon 
is of fundamental importance to the state and evolution of the 
universe, as outlined, for example, in the review by Heckman 
(1998). According to him, in the local universe starbursts are 
responsible for about 25% of the high-mass star formation rate 
and for 10% of the total luminosity; in the early universe, be- 
yond z ~ 0.7, IR luminous galaxies dominate the star forming 
activity (Flo'ch et al. 2005). Starbursts are also cosmologically 
significant if one interprets the high bolometric luminosities of 
high redshift galaxies to be due to star formation (Hirashita et 
al. 2003) at a rate so high that it can only be maintained over 
a cosmologically short spell (< 10** yr). The study of nearby 
starbursts would then help us to understand the processes un- 
derlying the star formation history of the universe. 
Starbursts, we think, are triggered by the gravitational interac- 
tion between galaxies (Kennicutt et al. 1987), but they occur, 
as a result of mass and angular momentum transfer, predomi- 

Send offprint requests to: rsiebemn@eso.org 



nantly in their nuclei, at the center of a massive and dynami- 
cally relaxed cluster of old stars (the bulge). Although the re- 
gion where OB stars form is relatively small (a few hundred 
parsec), its luminosity often exceeds that of the host galaxy. 
Starbursts are almost pure infrared objects, opaque to stellar 
photons. Whereas, on average, in the local universe ~60% of 
the star formation is obscured by dust (Takeuchi et al., 2006), 
in starbursts the fraction is typically 90%. To interpret infrared 
observations and to arrive at a self-consistent picture for the 
spatial distribution of stars and interstellar matter in the star- 
burst nucleus and of the range of dust temperatures, one has to 
simulate the transfer of continuum radiation in a dusty medium. 
Line emission is energetically negligible. 

A starburst has four basic parameters; total luminosity, L, dust 
or gas mass, Md or Mgas, visual extinction. Ay, and size. Size, 
Ay and Md are, of course, related, for a homogeneous density 
model, only two of them are independent. The luminosity fol- 
lows observationally in a straight forward way by integrating 
the spectral energy distribution over frequency, Md is readily 
derived from a millimeter continuum data point, if available, 
and the outcome is almost independent of the internal struc- 
ture or viewing angle of the starburst. The size is best obtained 
from radio observations as it does not suffer extinction by dust. 
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It does not matter whether the radio emission is thermal or non- 
thermal, both are connected to high-mass stars. 
In this paper, we present a set of SEDs for starbursts covering a 
wide range of parameters. Anyone with infrared data and inter- 
ested in their interpretation can compare them with our models, 
find an SED that matches (after normahzation of the distance) 
and thus constrain the properties of the starburst under investi- 
gation without having to perform a radiative transfer computa- 
tion himself. 

2. Dust model and radiative transfer 

A description of the dust model and the radiative transfer can 
be found in chapter 12 and 13 of Kriigel (2003), we here only 
summarize the salient points. We use standard dust. It consists 
of silicate and amorphous carbon grains with an MRN size dis- 
tribution (n(fl) oc fl"3^, a ~ 300 . . .2400A), and a population 
of small graphite grains (a ~ 10 . . . lOOA). There are also two 
kinds of PAHs (A^c = 38,A^h = 12 and A^c = 250, A^h = 48, 
where A^c^ ^'h are the number of C, H atoms, respectively). By 
mass, 63% of the dust is in siUcates, 37% in carbon of which 
60% is amorphous, 38% graphitic and 2% in PAHs. About 5% 
of the graphitic particles are so small (< 60A) that their tem- 
perature fluctuates. Such a dust mixture produces reddening in 
rough agreement with the standard interstellar extinction curve 
for/?v = 3.1. 

An important feature of our model is the division of the sources 
in the starburst nucleus into two classes. 

a) OB stars in dense clouds with total luminosity Lqb- The 
inmiediate surroundings of such a star constitutes a hot spot and 
its emission must be evaluated separately as they, more than 
anything else, determine the MIR part of the SED of a galac- 
tic nucleus (Kriigel & Tutokov 1978, Kriigel & Siebenmorgen 
1994). The outer radius of a hot spot, /?hs is given by the condi- 
tion of equal heating of the dust from the star and from the am- 
bient radiation field. The hot spots, whose total volume is small 
compared to the volume of the galactic nucleus, are presented 
in the radiative transfer equation by a continuously distributed 
source term Sy'^ir), where r is the distance towards the center 
of the galactic radius. For a fixed OB stellar luminosity, e^* is 
sensitive to the assumed density in the hot spot, 

b) All other stars of total luminosity Ltot - Lqb. These are 
mainly the old bulge stars of low brightness and surface tem- 
perature, but also hotter stars not enveloped in a dense cloud. 
This population is presented in the radiative transfer equation 
by a continuously distributed source term ev"'^'^(r). 

The model galactic nucleus is a sphere (of radius R) and the 
radiative transfer is computed with ray tracing. The intensity, 
^v(p. z), is a function of frequency v, impact parameter p, and 
coordinate z. At different v and p, we solve along the z-axis the 
equations 



rit) 



!' 

Jo 



/+(t) = /+(0) e"^ + r 5(x)e-'-^dx 
Jo 



Six)e''~'dx 



(1) 



(2) 



is zero at z = and increases with z, the optical depth t is 
zero at the edge of the nucleus (where Ze = -\lR~ - p^) and 
decreases with z. There is no radiation incident from outside, 
so I~{z = Ze) - l^{t - 0) = 0, and symmetry requires = I~ 
at z = 0. The source function (dropping sums over different 
kinds of dust particles) equals 



Sy = 



+ fiv"'^" + /sTf + K^"' I PiT)ByiT) dT 



(3) 



and / refer to the plus and minus direction of z, respectively. 
The indices p and v have been omitted. The optical depth r 



The term K'f'^ j P(T)By(T) describes the emission of dust 
grains. If they are big, the probability density PiT) equals the 
(5-function SiTa) where follows from the equilibrium be- 
tween radiative heating and cooling. For small grains, P(T) 
is evaluated in an iterative scheme similar to the method of 
Guhathakurta & Draine (1989). Jy is the galactic radiation field. 
As we assume isotropic scattering, we reduce the Mie scatter- 
ing efficiency, 2*"^", by the factor (1 -gv), where gy is the asym- 
metry factor. All quantities depend on the galactic radius, r. 
The emission from the hot spots is calculated separately in a 
radiative transfer program for an OB star in a spherical cloud 
of density p^^ bathed in the galactic radiation field. 

3. Parameter space of the model grid 

For our set of SEDs, we vary in the calculations the following 
five parameters: 

1. total luminosity L"" from 10'° to IO'^'Eq in steps of 0.1 in 
the exponent, 

2. nuclear radius, R = 0.35, 1 and 3 kpc, 

3. the visual extinction from the edge to the center of the nu- 
cleus: Av ^ 2.2, 4.5, 7, 9, 18, 35, 70 and 120 mag, 

4. ratio of the luminosity of OB stars with hot spots to the total 
luminosity: LoB/i'°' = 0.4,0.6 and 0.9, 

5. dust density in the hot spots. For a gas-to-dust ratio of 150, 
the corresponding hydrogen number densities are n^^ = 
10^, 10^ and lO'' cm"^. The density is constant within the 
hot spot. 

Not all parameter combinations are included in the set of SEDs 
because some are astronomically unlikely (for instance, very 
high L'°' and very little extinction). Altogether, the grid con- 
tains 7000 entries. 

The dust density in the nucleus, p, is spatially constant, dp/dr - 
0. Its value follows from the extinction Ay and the nuclear 
radius, R. The dust mass, Mj, is then given by AnpR^ IT> and 
increases linearly with Ay. For example, for R - 350pc and 
Ay = 18, the gas mass, Mgas, is 1.7 x IO^Mq. The density of 
all stars is centrally peaked, p«(r) oc r"'-^. The OB stars are 
always confined to the inner 350 pc, and they have a fixed lu- 
minosity and surface temperature (2x10* L©, T^s = 25000 K). 
The bulge stars fill the total volume. As they do not form hot 
spots, we need not specify the luminosity of a single star. For 
Ltot ^ IO'^-^Lq, their surface temperature, T^s, equals 4000 K 
(old giants). When Ltot > lO'^'^Lg, we assume T^s = 25000 K, 
which means that they consist mainly of OB stars, but outside 
compact clouds. 
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Fig. [2 illustrates the changes in the SED when one parameter 
is varied while all others stay fixed; fluxes refer to a source dis- 
tance of 50 Mpc. Panel a) informs us how a rise in luminosity 
shifts the far IR peak to shorter wavelengths. The flux then in- 
creases in the near IR much more strongly than at submillime- 
ter wavelengths. If Ltot > lO'^'^L©, the large grains become 
so warm that at A > I Ifim they outshine the PAH features. In 
panel b), we see that as the source becomes bigger, the dust 
gets cooler (maximum emission at longer wavelengths). This 
is not because the dust is then, on average, farther away from 
the source, but because the dust mass, Md, grows with when 
Ay is constant, and the mean dust temperature is determined 
by L'°7Md. We also see that a high ratio of Lqb/L}"^ enhances 
the near IR flux. Panel c) shows the influence of the density 
in the hot spots for n^^ - 10^, 10"^ and 10"^ cm"^, it is partic- 
ular strong in the MIR. Panel d) depicts the influence of the 
optical depth. Large values suppress the near IR emission and 
produce absorption in the lOyum, and for very high extinctions 
(Ay > 70mag) also in the IS/im siUcate bands. 

4. Testing the SED library 

4.1. Fitting prototypical galaxies 

We put our library of ~7000 theoretical SEDs to the test 
by applying it to seven famous and well studied galaxies of 
the local universe, five starbursts (M82, NGC253, NGC7714, 
NGC1808, NGC7552) and two ULIRGs (NGC6240 and 
Arp220). The observational data and our fits are displayed in 
Fig|5]to Fig|5] underlying model parameters are listed in Tab. 1 . 
The seven galaxies, discussed in more detail below, cover a 
wide range of luminosities and we check: 

- whether their observed spectra can be reasonably matched 
or, at least, bracketed by elements of the set; 

- whether there is only one matching element or, at most, a 
few which are similar in their basic parameters; 

- whether the parameters of the matching element are 
meaningful, i.e. whether they are consistent with the 
information about the structure of the galactic nucleus 
which we already have. 

M82 

The present model for this archetype starburst is similar to the 
one proposed before (Kriigel & Siebenmorgen 1994). The latter 
was shown only for A > 3jum. At shorter wavelengths, the ob- 
served flux does not steeply decline, as the old model predicts 
and as one would expect judging from the silicate feature (its 
depth implies Ay > 15 mag). Therefore, either hard radiation 
leaks out because of clumps or funnels created by supernova 
explosions, or there are stars in M82 outside the opaque nu- 
clear dust clouds. As our model cannot handle clumping, but 
we nevertheless wish to extend the spectrum into the UV, we 
simply add another stellar component. It is not included in a 
self-consistent way, but as its luminosity is ~ 10% of the total, 
such an approximation may be tolerable. The stellar tempera- 
ture and foreground reddening of the additional component are 



poorly constrained (see caption of Fig|2li. This is also reflected 
by the controversial interpretations via an old stellar population 
(Silva et al. 1998) or via young, but obscured stars (Efstathiou 
et al. 2000). 

We mention that Sturm et al. (2000) contest the existence of the 
lOjum silicate absorption feature in M82. They think that the 
value T(18jum) / T(9.7yum) is too low and that recombination 
line ratios, like H^/Ha, indicate only Ay ^ 5. However, the 
small ratio T(18;um)/T(9.7yum) is a radiative transfer effect 
where 18/im emission is favored over lOyum emission, and 
the H^/Hq, ratio increases when the dust is not a foreground 
screen, as Sturm et al. assume, but mixed with the HII gas. The 
strongest argument that the lOyum depression is due to silicates 
comes from the 1mm flux. It implies a large column density 
of cold dust which, unless it is all behind the hot dust, must 
produce the absorption feature. 

NGC253 

NGC253, another bright and nearby starburst, shows at high 
resolution in the MIR a complex structure with several knots 
(Galliano et al. 2005). Nevertheless, the low spatial resolution 
observations are well reproduced in our fit (Fig|3l which, in 
the 10 - 40jum wavelength range, is of similar quality as in 
M82 (FigE). Below 2 Jy, ISOSWS data are noisy and have 
therefore been omitted. The dip at 18;um in the model of 
Piovan et al. (2006) is not present in ours, and not borne out 
by the observations. As our model dust is uncoated, the ice 
features reported by Imanishi et al. (2003) are not reproduced. 

NGC7714 

Spitzer spectra (Brandl et al. 2004) do not reveal AGN 
signatures and support our interpretation that the nucleus is 
dominated by a weakly obscured starburst (FigO. 

NGC1808 

This starburst is claimed to be young (Krabbe et al. 1994). In 
high resolution MIR images, several hot spots are detected 
which coincide with the most intense radio sources (Galliano 
et al. 2005). Our fit (Fig|3 somewhat underestimates the PAH 
emission. So one may have to increase the PAH abundance (as 
was done in the model of Siebenmorgen et al. 2001) which is 
set constant in the computations of our SED library. Piovan et 
al. (2006) predict silicate absorption features at 10 and 18yum 
which, however, are not detected. The dip at ~ lOfim is due 
to the wings of neighboring PAH bands. It is not caused by 
silicate self-absorption which would require much higher 
optical depths. 

NGC7552 

This infrared luminous galaxy harbors a ring-like circumnu- 
clear starburst (Siebenmorgen et al. 2004). Neglecting such 
structural details, our fit to the dust emission is satisfactory. 
Two models are shown in Fig. |3 which bracket available data. 
The hot spot density is low and the OB luminosity ratio, 
LoB/Ltot, is not well constrained (Tab. 1). 

NGC6240 
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Table 1. Fit parameters to models of Figs. 2-5 
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NGC6240 is a merging ULIRG. Such objects are one to two 
orders of magnitude brighter than starbursts. Our fit in Figls] 
is acceptable despite a ~ 30% deficiency near AQ-SQjjm. To 
better match the NIR photometry, we added to the starburst 
SED a 4000K black body with L = lO'^ '^Lo. Lutz et al. (2003) 
suggest that stars account for most (-75%) of the total lumi- 
nosity and that the rest is due to an optically obscured AGN. 
When they subtract from the SED of NGC6240 a scaled-up 
M82 template, a faint (0.07Jy) residue remains which they 
attribute to the AGN. Dopita et al. (2005) underestimate in 
their model the 10-30yum region (by a factor ~4 at 15jt/m) 
but they argue that they could match the data if they added 
an AGN component, a procedure which is sometimes applied 
to galaxies with hidden broad line regions (Efstathiou & 
Siebenmorgen 2005). 

Arp220 

Arp220 is the nearest example of a ULIRG. MIR high reso- 
lution maps (Soifer et. al 2002) show a double nucleus with 
1" (360pc) separation. We process low resolution Spitzer IRS 
data using the SST pipeline (Higdon et al. 2004). The ISOPHT 
(Spoon et al. 2004) and Spitzer spectrum reveals a complex 
spectrum with ice and silicate absorption and pronounced R4H 
emission bands at 6.2 and 7.7/im. Dopita's et al. (2005) model 
predicts PAH features that are too strong (factor > 5). Piovan et 
al. (2006) fit the SED of the central 2kpc region using an optical 
depth of Tv = 35mag and a dust model with an SMC extinc- 
tion curve. Siebenmorgen et al. (1999) proposed ry - 54mag 
and MW dust to fit the photometric data available at that time. 
The SED library fit gives R - 3kpc and Ay =72mag. A model 
SED with Ay =120mag and R - Ikpc yields too strong silicate 
absorption and requires an additional cold dust component for 
the submm. 

4.2. Predicting fluxes 

The infrared luminosity is the key parameter of a galaxy and 
it is often used to estimate the star formation rate of a galaxy 
(Kennicut 1998). Unfortunately, for faint or redshifted objects 



Fig. 6. All elements of the SED library (dotted) which fit 8 and 
24jum photometry (circles) of NGC6240 to within 30%. Best 
fit (full line) and other data as of Fig|5l 

photometry is sometimes only provided at two MIR bands, for 
example at 8 and 24yum from the Spitzer satellite. 
In Fig. |6lwe demonstrate for the ULIRG NGC6240 that the 
SED library can be used to estimate the total luminosity to 
within a factor ~ 2 from two such MIR fluxes only. We also 
see (Fig.|6} that the SED will be quite well constrained by an 
additional submm data point. 

5. Discussion 

5.1. Methods of modeling starburst SEDs 

One finds in the literature three different ways to reproduce or 
explain the SED of an extragalactic object. 

i) Matching it with a template SED of a well known galaxy 
(Laurent et al 2000, Lutz et al. 2003, Spoon et al. 2005). This 
is reasonable only as long as template and object are similar 
in their parameters as well as geometrical structure and orien- 
tation on the sky: So it is not meaningful to compare SEDs 
of AGN type 1 and type 2, or objects with radically different 
luminosities, like M82 and NGC6240, because the luminosity 
aff'ects the SED, as can be seen from Fig.[^. To obtain a good 
match between the SED of the object and the template, after 
normalization to a unit distance, one usually has to scale the 
flux of the template moderately up or down to fine-tune the lu- 
minosity. If the fit is successful one may argue that the object 
is similar to the template in its geometry and basic parameters 
(but for some scale factor) and that one understands it almost 
as well as the template. 

ii) Reproducing the shape of the SED by optically thin dust 
emission. The dust is assumed to be heated in a given radiation 
field which is usually a scaled-up version of the interstellar 
field (Devriendt et al. 1999, Dale et al. 2001, 2005, Lagache et 
al. 2003). This procedure neglects all effects of radiative trans- 
fer and must fail when dust self-absorption becomes important, 
most strikingly in the ~ lOjum region as shown in Fig.^. 
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Hi) Solving the radiative transfer, in various degree of so- 
phistication. This is a much more ambitious method and re- 
quires assumptions about the structure and parameters of the 
galaxy. Three-dimensional codes have been applied to spi- 
rals (Kylafis & Bahcall (1987), Popescu et al. (2000), Tuffs 
et al. (2004)) using ray tracing or Monte Carlo techniques 
(Bianchi et al. 2000). Rowan-Robinson & Crawford (1989) fit 
IRAS color diagrams of starbursts using a one-dimensional 
transfer code. 

We also do radiative transfer calculations and it may be instruc- 
tive to point out technical and conceptual differences between 
our models and those devised recently by other authors (Silva et 
al. (1998), Efstathiou et al. (2000), Takagi et al. (2003), Dopita 
et al. (2005) and Piovan et al. (2006)), although we admit that 
we did not always find it easy to pin down exactly which ap- 
proximations our colleagues used (as they may experience dif- 
ficulties in identifying our assumptions). 
All groups evaluate the emission from a dusty medium of 
spheroidal shape filled with stars, all seem to use similar optical 
dust constants, and all incorporate small grains with tempera- 
ture fluctuations (like PAHs). At first glance, the model results 
appear to agree, but upon closer inspection one finds that the 
maximum deviations in the other papers are considerable (fac- 
tor 4), whereas the fits from our library to the prototype objects 
(Fig. 2 to 5) are nowhere off by more than 50%. 

The major points where our paper differs concern the treatment 
of the sources, the interstellar extinction curve, the radiative 
transfer, and the presence or neglect of hot spots. 

a) Stellar sources. We do not take into account the time 
evolution of a stellar population after the burst and the possi- 
bility that there may have been several episodes of rapid star 
formation. Our models are therefore simpler and do not allow 
to constrain the age of the burst(s). We are, however, skepti- 
cal that this is possible without spectroscopic observations (like 
stellar CO bands in supergiants), also derived ages of the stellar 
populations are controversial (see M82). 

b) Extinction curve. We assume galactic dust and do not 
consider the possibility that it may be a combination of the 
species found in the Milky Way, LMC or SMC. Again, here 
our model is simpler, but as the extinction of the sources is 
usually large (Ay > 5 mag), the exact shape of the reddening 
curve has little effect on the resulting infrared SED. 

c) Radiative transfer As far as we can tell, the interac- 
tion between dust and radiation is treated consistently only in 
Takagi et al. (2003) and the present paper Other authors intro- 
duce basic emission units which they compute separately. Such 
a unit may be a diffuse gas cloud, or a spherical molecular cloud 
filled with dust and stars that are either continuously smeared 
out over the cloud or concentrated in the center. The emission 
units are then scaled up, or a simplified radiative transfer (with 
constant source function or no reemission) is applied to match 
the nucleus under consideration. Naturally, when the optical 
thickness is not small, models without radiative transfer are at 
some point faulty, although it is hard to quantify how much the 
simplifications effect the resulting SED. 

d) Hot spots. They are a particular feature of our mod- 
els and inevitably arise when a luminous star is enveloped by 



a cloud with a density considerably above the mean density 
of the nucleus. Neglecting hot spots seriously underestimates 
the MIR emission of the nucleus (see Kriigel & Siebenmorgen 
(1994)orFig.[Il;). 

e) dumpiness. As discussed in the model description of 
M82, the optical and UV flux is best explained by postulating 
the interstellar medium to be clumped. Clumping is a natural 
consequence of SN explosions; when the surface filling factor 
is close to one, it has little effect on the SED at wavelengths 
greater than a few micron. The models that include stellar evo- 
lution (Silva et al., Efstathiou et al., Takagi et al., Dopita et 
al., Piovan et al.) introduce as an additional free parameter the 
fraction of starlight escaping the galaxy due to clumping; this 
fraction depends in their computations on the age of the star- 
burst. Our approach is again simpler. Because the UV and op- 
tical stellar light that leaks out is in reality modified in a com- 
plicated way by the passage through a clumped medium, we 
only add, where necessary, a blackbody curve to account for 
the excess light. 

5.2. Completeness, uniqueness and credibility of the 
SED library 

It is remarkable that one can very well fit the SEDs of galaxies, 
like M82, Arp220 and others, with models of constant density 
and radial symmetry. The satisfactory fits imply, first, that our 
library grid is sufficiently fine, and we expect that starbursts 
observed with similar wavelength coverage as those presented 
in Fig. 2-5 (more than a few data points in the SED) can be 
reasonably matched by a single element of the SED library. 
Nevertheless, one may wonder whether the fits are meaning- 
ful. After all, we know for AGN, which have tori that lead to 
the division into type 1 and 2 with respect to the observer, that 
spherical symmetry is a principally unacceptable approxima- 
tion. As the torus is the result of rotation, it should form inde- 
pendently of a massive black hole and therefore also exist in 
starbursts. However, there seems to be no need to invoke one. 
There are probably two explanations. First, whereas an AGN is 
small (pc) and easily shadowed by the much bigger torus (100 
pc), a starburst region is as large as or larger than a torus and 
could not be blocked visually. So there cannot be starbursts of 
type 1 and 2. Second, the galaxy collision preceding the star- 
burst leads to strong perturbations of the nuclear gas which, 
in the gravitational potential of the little disturbed bulge stars, 
results in rough spherical symmetry. 

We also have to discuss the possible contamination of starburst 
fluxes by emission from the galactic disk when the spatial reso- 
lution of the observations is poor. This is the standard situation 
in the far infrared. There, however, the nucleus is usually much 
brighter than the disk and then the contamination is irrelevant. 
It may be substantial at short wavelengths (NIR, optical, UV) 
if the starburst is very obscured and little optical flux leaks out. 
In the submm/mm region, one measures mainly the dust mass 
and there is likely to be more mass in the disk than in the 
core. To estimate the contributions of the disk and the core, let 
F, L, M and T denote the observed flux, bolometric luminosity, 
dust mass and dust temperature, respectively. With the approx- 
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imations L oc MT^ and Fimm °^ MBx^^iJ) oc MT, the flux ra- 
tio, say, at 1 mm, becomes f imm^c/f imm.d = {l-cll-dWilT^f . 
Here we used the index d for disk and c for the core. Typical 
mean values are - 10...20K (Kriigel et al. 1998) and 
Jc = 30 . . . 50 K (Klaas et al. 2001). Therefore, the cold dust in 
the disk is not really important in measurements of low spatial 
resolution as long as the nucleus is much brighter than the disk. 

6. Conclusions 

We have computed in a self-consistent radiative transfer SEDs 
of spherical, dusty galactic nuclei over a wide range of their 
basic parameters such as luminosity, dust mass, size and ob- 
scuration. The SEDs can be accessed in a public library ' . 
Given a set of data points for a particular galaxy, there is a sim- 
ple procedure, described in the README file, to select from 
the library those elements which best match them. If the obser- 
vations cover the full wavelengths band from a few to about 
1 mm, one usually finds only one library element that fits very 
well, as demonstrated for seven famous active galaxies. If the 
data points are widely spaced, there may be a few elements of 
less fitting quality, but similar in their basic parameters. 
The library therefore allows one to constrain the fundamental 
properties of any nucleus which is powered by star formation 
and for which data exist, without any further modeling. Two 
observed fluxes in the MIR plus one submm point are usually 
sufficient for a crude characterization of the nucleus. If there 
are only two MIR points, for example, from Spitzer at 8 and 
24/im, one can still bracket the total luminosity within a factor 
of ~2. 

In the UV, optical and NIR, it may be necessary to add to the 
library SED a low luminosity stellar component, at least, this 
was necessary for M82 and NGC6240. This component may 
be due to photons that escaped the nucleus without interaction 
because of clumping or it may be light from the galactic disk. 
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Fig. 1. Influence of starburst parameters on the SED for a distance of 50 Mpc. The parameters which are kept constant are listed 
in square parentheses, a) Total luminosity is var ied between L'°' = 10^° and IO'^'^Lq; [R = 3kpc, Ay ~ 17mag, Lob/^'"' = 0.6 
and n^' = 10^^ cm ]. b) Here we vary two parameters: the radius of the nucleus from R - 0.35 over 1 to 3kpc, and the luminosity 
ratio: Lob/^'"' = 0.4 (fuU Unes), 0.9 (dotted); [L™ = IO'^Lq, Ay ~ 4.5mag and n^' = 10"* cm-"*], c) Variation of the hot 
spot density: n^' = 10^ 10^ and 10"* cm^^; [L*°' = W'-^Lq, R = 3kpc, Ay ~ 9mag, Lob/^'"' = 0.9]. d) Variation of the dust 
extinction: Ay -2.2, 4.5, 6.7, 9, 18, 35, 70 and 125mag; [L'°' = IO'^^Lq, R = 3kpc, Lob/^'"' = 0.9 and n*"^ = lO'^cm-^]. 
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Fig. 2. SED of the central region of M82, data points with Icr error bar. Full line: library model with parameters in Tab. 1 . To fit 
the data below 5fim, we added to the SED library spectrum a blackbody, either unreddened (T - 2500 K, full line), or reddened 
(T = 8000K, Ay - 4mag, dashed, or T = 5000K, Ay = 3mag, dotted). Full 0.4 - ISOOyum wavelength range (top), a zoom into 
the 12 - 34//m region (bottom). Data references (1300yum: Krugel et al. (1990); 1 100 and SOOjum: Hughes et al. (1990); 400//m 
Jaffee et al. (1984); FIR: Telesco & Harper (1980), Rieke & Low (1972), Rieke et al. (1980), Telesco & Gezari (1992); IRAS 
NIR photometry in 40" - 100" aperture: Kleinmann & Low (1970), Jarret et al. (2003), Aaronson (1977) and Johnson (1966) 
between 2.3-40.4jum ISOSWS specti'um: Sloan et al. (2003)). 
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Fig. 3. SEDs of NGC253 and NGC7714, data points with Icr error bar. Models: full lines, model parameters in Tab. 1. Data 
for NGC253 (ISOSWS: Sloan et al. (2003); ISOLWS and ISOPHT: Radovich et al. (2001); NIR: Rieke & Low (1975); IRAS; 
submm: Rieke et al. (1973), Hildebrand et al. (1977), Chini et al. (1984)). Data for NGC7714 (NIR: Spinoglio et al. (1995), Jarret 
et al. (2003); Spitzer IRS: Brandl et al. (2004); IRAS; ISOPHT: Krugel et al. (1998); 850/zm: Dune et al. (2000); 1.3mm: Krugel 
et al. (1998)). 
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Fig. 4. SEDs of NGC1808 andNGC7552, data points with Icr error bar. Models: full and dashed lines, model parameters in Tab. 1. 
Data for NGC1808 (NIR: Glass (1976), Jarretet al. (2003);1RAS; ISOPHT 160jum andlSOCAM spectroscopy: Siebenmorgen et 
al. (2001). Data for NGC7552 (MR: Glass (1976), Jarret et al. (2003); ISOCAM: Roussel et al. (2001); T1MM12: Siebenmorgen 
et al. (2004); Spitzer IRS of nucleus: Kennicutt et al. (2003); IRAS; submm: Stickel et al. (2004), Hildebrand et al. (1977)). 
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Fig. 5. SEDs of NGC6240 and Arp220, data points with Icr error bar Models: full and dashed lines, model parameters in 
Tab. 1. To match JHK photometry of NGC6240, we added a 4000K black body to the starburst model. Data for NGC6240 
(NIR: SpinogHo et al. (1995); ISOPHT and submm: Klaas et al. (1997, 2001); 350;um: Benford (1999); ISOCAM specti'oscopy: 
(Laurent et al. 2000); ISOPHT and ISOSWS: Lutz et al. (2003)). Data for Arp220 (2MASS: Jarret et al. (2003); IRAS; ISOPHT: 
Klaas et al. (2001); ISOCAM: Siebenmorgen & Efstathiou (2001); submm: Benford (1999), Rigopoulou (1996), Dunne et al. 
(2000), Carico et al. (1992), Chini et al. (1986); ISOLWS archive spectrum is scaled to match the ISOPHT photometry and 
Spitzer IRS spectrum (this work)). 



